Thermal Stability and Electrical Control of Magnetization of Heusler /Oxide Interface 

and Non-collinear Spin Transport of Its Junction 
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Towards next-generation spintronics devices, such as computer memories and logic chips, it is nec- 
essary to satisfy high thermal stability, low-power consumption and high spin-polarization simulta- 
neously. Here, from first-principles, we investigate thermal stability (both structure and magnetiza- 
tion) and the electric field control of magnetic anisotropy on Co2FeAl (CFA)/MgO. A phase diagram 
of structural thermal stability of the CFA/MgO interface is illustrated. An interfacial perpendicular- 
anisotropy, coming from the Fe-0 orbital hybridization, provides high magnetic thermal stability 
and a low stray field. We find an electric-field-induced giant modification of such perpendicular- 
anisotropy via a great magnetoelectric effect (the anisotropy energy coefficient /3~10 — 7 erg/V cm). 
Our spin electronic-structure and non-collinear transport calculations indicate high spin-polarized 
interfacial states and good magnetoresistance properties of CFA/MgO/CFA perpendicular magnetic 
tunnel junctions. 

PACS numbers: 85.75.Dd, 75.30.Gw, 75.80.+q, 31.15.Ar 
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Perpendicular magnetic tunnel junctions (p-MTJs), 
are leading spintronics devices for next-generation com- 
puter memories (STT-MRAM) and logic chips for their 
high thermal stability and good magnetoresistance (MR) 
properties However, traditional hard ferromagnetic 
(FM) materials, such as Ll -FePt, are not good elec- 
trode materials. It is because of the poor lattice matching 
with the MgO barrier and low spin-polarization, which 
would sharply reduce the MR ratio. Furthermore, con- 
versional hard ferromagnets are also inclined to exhibit a 
large damping constant(a), e.g., cv( L i o -FePt)>0.055 due to 
the strong spin-orbit interaction (SOI), and hence a high 
switching current density, i.e., high-power consumption. 
Experimentally, Ikeda, et al. proposed the perpendicu- 
larly magnetized CoFe (a soft ferromagnet) in ultra-thin 
CoFe/MgO/CoFe MTJs to be effective for achieving high 
thermal stability and high MR ratio (124%) 0. How- 
ever, the damping constant of the CoFe film increases 
sharply with the decreasing thickness. Recently, an 
electric-field-assisted magnetic switching is demonstrated 
in CoFe/MgO/CoFe MTJs, which can efficiently reduce 
the switching current density Q. It provides a more 
energy-efficient route and is compatible with current 
ubiquitous electrically-controlled semiconductor devices. 
Actually, electric control of magnetization has been 
widely studied (both experimentally and theoretically) 
in ferromagnetic thin filmsQ, FM/metal[3-@, FM /oxide 
interfaces, such as FePt/MgO@, 13, Fe/MgO[iOl and 
CoFe/MgO@, H3 ]. Although the coefficient of magnetic 
anisotropy energy, /?, (A M ca = P x E F , where A M ca 
is the change in the magnetic anisotropy and E F is the 



electric field.) of FM/oxide is one order of magnitude 
higher than that of pure Fe surfaces Q and FM/metal 
interfaces but it is still quite low (salO -8 erg/V cm) 
for practical low-power applications. That is to say, the 
electrically-controlled energy-efficiency is quite low now. 
Therefore, To find a FM/oxide interface with a giant 
modification of magnetization by an electric field is a 
matter of great urgency. Meanwhile, it is still absent of 
FM/oxide/FM junctions to satisfy high thermal stability, 
low switching current density (i.e., low-power consump- 
tion) and good spin transport properties all at the same 
time. 

In this Letter, we predict a great magnetoelectric effect 
at the interface of Co 2 FeAl/MgO (001). The anisotropy 
energy coefficient, /3, is one order of magnitude higher 
than that of traditional FM/MgO interfaces. Further- 
more, high thermal stability (both structure and magne- 
tization) is demonstrated in the FeAl|0 interfacial struc- 
ture. The great magnetoelectric effect (CFA/MgO) com- 
bined with a small Gilbert damping constant (Co2FeAl) 
and low stray field (perpendicular anisotropy) means a 
low critical current density in CFA/MgO-based MTJs. 
On the basis of this interfacial structure, we build a 
junction (Co2FeAl/MgO/Co2FeAl) which shows a high 
spin-polarized interfacial states, and our non-collinear 
spin-transport calculations show this p-MTJ has better 
MR properties (MR ratio and ARA) than conventional 
CoFe/MgO/CoFe and FcPt/MgO/FcPt MTJs. The high 
thermal stability, low-power composition, and good MR 
properties indicate CFA/MgO/CFA MTJs to be one of 
the best candidates for next-generation low-power/high- 
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FIG. 1. (Color online) (a) Schematic structure and high-resolution cross-section transmission electron microscopy images of 
Co2FeAl/MgO(001) [ref.25]. (b) and (c) are interfacial structures of CFA and MgO with an indication of the size of unitcell 
(the square of white dash line). The different contacted sites are labeled by number 1-4. 



performance spintronics devices. 

A CFA/MgO(001) interface was modeled as shown 
in Fig. la. The in-plane lattice parameter was fixed 
at l/y/2 of the experimental lattice parameter of bulk 
CFA (5.730 A) with a very small lattice mismatching 
with MgO (~3.7%). The electronic structure calcula- 
tions of CFA/MgO (001) interfaces were carried out using 
a first-principles approach[scc Supplemental Materi- 
als]. The non-equilibrium Green's function formulism 
was used for the non-collinear transport calculation of 
CFA/MgO /CFA p-MTJs [see Supplemental Materi- 
als]. 

Recently, a half-metallic Heusler compound (Co2FeAl) 
is demonstrated the perpendicular magnetic anisotropy 
(PMA) in CFA/Pt multilayers prepared on MgO (001) 
substrates (l3j|. Subsequently, Wen el al. reported a 
PMA in the structure of CFA/MgO with a high K u 
of 2-3xl0 6 erg/cm 3 [l4|. Before any magnetoelectric 
or transport studies on CFA/MgO structures, the first 
and important thing is to determine the interfacial ther- 
mal stability structurally and magnetically. For exam- 
ple, although the high-resolution transmission electron 
microscopy (HRTEM) can give a clear atomic structure 
of the CFA/MgO interface, it is still not clear the dis- 
tribution of different types of atoms at the interface (see 
Fig. la). Furthermore, it is known that the interfa- 
cial geometry strongly affects the spin interfacial states 
and spin transport properties 0, [IH [HI], while only a 
few studies are carried on the the interfacial structure 
of ternary Heusler alloys with oxides because it is much 
more complex than that of CoFcB after high-temperature 
annealing [17Hl9l ]. Therefore, we first study the ther- 
mal stability of the interfacial geometry of CFA/MgO. 
Figure lb and lc show that the CFA/MgO interface 
has four possible types of CFA-terminations, i.e., CoCo, 
FeAl, FeFe, and A1A1. Furthermore, each of these termi- 
nations has four possible positions relative to MgO (001), 
i.e., O-top, Mg-top, hollow and bridge site. The total 
energy calculations show that for the CoCo-, FeAl-, and 
FeFe-terminated CFA, the configuration of the O-top site 



is lower in energy than any others, whereas for the A1A1- 
termination, the Mg-top site is the most energetically fa- 
vorable. Thus, we use the method of ab initio atomistic 
thermodynamics to compare the stability among differ- 
ent terminations 17 , 2^] . The interface formation energy 



is calculated by the formula: 
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where AG is the difference in the Gibbs free energy 
between the heterostructure and appropriate amounts of 
two bulk materials forming the interface. ANi (i = Co, 
Fe or Al) counts the number of atoms, in which the su- 
percell deviates from bulk composition, and fa represents 
the chemical potential of the element i. Since the inter- 
face energy is a free-energy difference, it is often a good 
approximation to calculate 7 from the differences of to- 
tal energies obtained from DFT calculations [2l|. For the 
thermodynamic analysis, it is assumed that the inter- 
face is in thermodynamic equilibrium with bulk Heusler 
compounds. This allows us to eliminate one chemical po- 
tential term with the relation: c/cfa = 2/i,co + MFc + Mai, 
where <7cfa represents the Gibbs free energy of Co2FeAl 
per formula unit. Hereby, only two of the three chemi- 
cal potential values fa are independent. Under the Co- 
or Fe-rich condition, Co or Fe is assumed in thermody- 
namic equilibrium with its bulk solid phase (/ico(buik) 
and /LtFe(buik) ) because the excess Co or Fe atoms may 
lead clusters or even precipitation. Consequently, the 
chemical potential uq and /ipc may not exceed their 
chemical potential of bulk phase, uqo < /UCo(buik) an d 
MFc < MFo(buik)- The interface phase diagram (Fig. 2) 
displays the termination with the minimal interface en- 
ergy for every chemical potential pair (fico and MFc)- It 
can be concluded from the phase diagram that two ter- 
mination configurations of the Heusler compound, i.e., 
the CoCo- and FeAl-planes sitting on top of oxygen, are 
thermodynamically stable. In contrast, the other two 
terminations are not accessible within the limit set by 
thermodynamic equilibrium with other phases. Further 
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FIG. 2. (Color online) Phase diagram for the epitaxial 
CFA/MgO (001) interface. The colored regions correspond 
to different interface terminations being stable under the con- 
ditions described by the chemical potentials A^i; = fit — gt, 
where (i =Co,Fe) is the total energy of one unitcell. The 
orange hex-polygon indicates the region accessible in thermo- 
dynamic equilibrium with the bulk phases of Fe, Co, FeaSi, 
Co 3 Si, FeSi, or CoSi. 



analysis on the phase diagram shows that, in the wide 
range of ^co, the FeAl-termination is more favorable than 
the CoCo-termination because almost all the orange area 
are contributed by the green area. In experiments, the 
C02FCAI film is high-temperature annealed after deposi- 
tion on a buffer layer in order to obtain the ordered L2i 
phase. Therefore, the most thcrmodynamically stable in- 
terface structure after annealing in the experiment should 
be the one derived from our chemical potential analy- 
sis above, i.e., FeAl|0. Thus, we following focus on the 
FeAl-terminated Co 2 FeAl/MgO(001) interface and dis- 
cuss its magnetocrystalline anisotropy, magnetoelectric 
effect and spin transport. 

Besides the structurally thermal stability, the magnet- 
ically thermal stability (hard magnetism with PMA) at 
the interface plays another key role on the device per- 
formance because the perpendicular magnetic anisotropy 
{K u ) can beat the thermal perturbation in ambient con- 
ditions (k-QT, where, fee and T arc Boltzman constant 
and temperature, respectively). Without PMA, the MR 
ratio of MTJs in room temperature would be much 
lower than that in low temperature due to the ther- 
mal perturbation on the interfacial magnetization pjj . 
Moreover, the perpendicular anisotropy can strongly re- 
duce the stray field, and then the power consumption 
of devices. For CFA/MgO, our calculations (involv- 
ing SOI) show that perpendicular magnetocrystalline 
anisotropy is energetically favored with an anisotropy 
constant K u of 0.428 erg/cm 2 , which is in good agree- 
ment with the previous experimental report |14j. The ori- 
gin of PMA can be understood by the electronic structure 
around the Fermi level via the second order perturbation 
consideration 0, [H, 23 1: 
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FIG. 3. (Color online) (a) The orbital-resolved PDOSs pro- 
jected on the interfacial iron and oxygen atoms in absence 
of electric fields (dash line) and under the electric field of 
2 V/nm (solid line). Induced differential charge density, in 
units of e/A 3 , at the interfacial Fe atom for E=2 V/nm in 
the (100) plane (b) and (001) plane (c). 
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where MCA is the magnetocrystalline anisotropy; k 
and k, u represent the occupied and unoccupied states 
with the wave vector k; and L z and L x are the angu- 
lar momentum operators along (001) and (100) direc- 
tion. According to Eq. 2, the SOI between the oc- 
cupied and unoccupied states with the same magnetic 
quantum number (771) through the L z operator enhances 
£-MCA(k), while that with different m through the L x op- 
erator weakens it. Figure 3 shows the projected density 
of states (PDOS) of CFA/MgO. Only the minority-states 
are shown for clarity, since the MCA originates mainly 
from the SOI between the minority spin bands around 
the Fermi level 0, H, [23|. A metal-induced gap state 
(MIGS) is occurred near the Fermi energy, which is com- 
posed of 0-p z and Fe-e2 g orbitals (d z 2 and d x i_ y i). Thus, 
on the basis of the above perturbation theory, the spin- 
orbit coupling between 0-p z and Fe-d z 2 (with the same 
quantum number 0) through the L z operator (001) con- 
tributes positively to the PMA in sense that the matrix 
element < p z \L z \d z 2 > is non- vanishing, and hence ex- 
plains the perpendicular (001) magnetic anisotropy. Note 
that the contribution to MCA of < p z \L z \d x 2_ y 2 > and 
< Pz\L x \d x 2_ y 2 > are zero. 

Recently, a manner of efficiently manipulating vari- 
ous magnetic properties by the application of an ex- 
ternal electric field, namely a magnetoelectric effect. 
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FIG. 4. (Color online) The magnetic anisotropy energy den- 
sity (t=l nm) of CFA/MgO as a function of the external 
electric field. Some experimental and calculated K u ■ t of 
CFA/MgO, CoFe/MgO and Fe/MgO (without E F ) are shown 
for comparison. The inset table shows the anisotropy coeffi- 
cient (/3) of some FM/MgO structures (with E F ). 
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FIG. 5. (Color online) The spin-resolved LDOSs at the Fermi 
energy projected on to each atomic sphere as a function of 
the distance from the interface for Co2FeAl/MgO/Co2FeAl 
(001). The positive (negative) sign of y axis indicates spin up 
(down) . 



has been successfully demonstrated theoretically and 
experimentally P, 0-0 Q. As mentioned in the introduc- 
tion, electrical control of magnetization is not only more 
power saving, but also compatible with current ubiqui- 
tous electrically-controlled semiconductor devices for the 
large-scale integrated circuit design. Thus, we next inves- 
tigate the electric-field-induced modification of MCA of 
CFA/MgO. In Fig. 4, our calculated anisotropy energy 
density of CFA(i =1 nm) /MgO is in good agreement with 
previous experimental results [3] and in the same order of 
magnitude with calculated CoFe/MgO[l| or Fe/MgOQ. 
Our calculated values of the MCA density show a almost 
linear dependence on the external electric field. The in- 
terfacial MCA coefficient (3, defined as AK U ■ t = /3 X E F , 
where AK is the change in the MCA energy, a key index 
for the strength of the MCA modification by an electric 
field. Surprisedly, the value of (3 of CFA/MgO is one or- 
der of magnitude larger than that of CoFe/MgO[l],[l2j or 
Fc/MgO[l[Iol and almost two order of magnitude larger 
than that of Fe thin film or Fc/metal (001)0, i, [2J]. 
Such giant modification of MCA by the electric field in 
CFA/MgO indicates a more energy-efficient route from 
the electric control of magnetization for low-power de- 
vice applications. 

It is found that the above great ME effect derives 
from the electrical-ficld-induced redistribution of spin 
electrons between difference orbitals around the Fermi 
level. In Fig. 3a, the originally degenerated high-lying 
Fe-e2 9 states are split by the external electrical field with 
an exchange splitting energy of 0.53 eV. Such splitting- 
pushes the ¥e-d x 2_ y 2 orbital below the Fermi level of 
0.4 eV (thus fully occupied). Meanwhile, the 0-p z or- 
bital follows up the shift of Fe-d x 2_ y 2 orbital and always 



strongly couples with Fe-d x 2_ y 2 independent on the ex- 
ternal electric field. The Fe-d z 2 is shifted above the Fermi 
level by the electric field, indicating a charge transfer 
from Fe-d z 2 to Fe-d x 2_ y 2 and 0-p z . Consequently, the 
mentioned origin of PMA (< p z \L z \d z 2 >) is completely 
destroyed due to the decoupling between 0-p z and Fc-d z 2 
[see Fig. 3a]. A further quantitative illustration of such 
d-electron redistribution is shown in Fig. 3b and 3c by 
plotting the induced charge density (differential) on in- 
terfacial Fe atoms, ApFc = PFc(2.0) — pFc(O), within the 
energy window Ef— 0.5 eV < E < Ei + 0.2 eV, where the 
specified charge transfer takes place. Figure 3b and 3c 
clearly show a reduced occupation of the Fe-d z 2 orbitals 
and an enhanced occupation of Fc-d x 2_ y 2 orbitals. 

We note that in the experiment the Li2-FePt/MgO 
interfacial structure also shows good interfacial proper- 
ties, such as large PMA (if„~10 6 -10 7 erg/cm 3 ) and gi- 
ant anisotropy energy coefficient ((3^10~ 7 erg/V cm)[7[. 
However, FePt/MgO/FePt p-MTJs show worse MR 
properties, neither TMR ratio nor AMR, because of low 
spin-polarized interface states [3]. This raises the ques- 
tion of whether CFA/MgO (with good interfacial prop- 
erties) has good transport (i.e., MR) properties for prac- 
tical device applications? With this question, we next 
investigate the spin-polarization of CFA/MgO/CFA p- 
MTJs, including interfacial spin polarization and non- 
collinear spin transport properties. 

It is well known that the spin-dependent tunneling 
conductance of the MgO-based MTJs is determined 
by two crucial factors. One is the symmetry match- 
ing between the transmission Bloch states in FM elec- 
trodes and evanescent states in the MgO barrier. The 
other is the spin-polarization of the FM/MgO interfacial 
states [l5, 3 3- The Bloch states with Al symmetry 
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TABLE I. The calculated conductance (Siemens) of the 
parallel majority (G™ a jL), parallel minority (G™™) and 
antiparallel {G™%/ min ) channels of the CFA/MgO/CFA, 
CoFe/MgO/CoFe, and FePt/MgO/FePt perpendicular mag- 
netic tunnel junctions. The MgO spacer length is 17.1 A for 
all three structures. 
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suffer from the smallest decay in MgO[l9j. Therefore, 
the spin-filter efficiency can be tremendously enhanced if 
the FM electrodes have Al-band half-metallicity. It has 
been reported that, for both bulk C FA (2^1 and CoFeflU 
electrodes, only majority bands which cut the Fermi 
level possess the Al feature and fulfill the symmetry- 
matching criterion. In contrast, Al-band of Ll2-FcPt 
predominantly contributes to the tunneling conductance 
for both the majority and the minority-spin channels. p6j 
Figure 5 illustrates the spin-resolved layered LDOSs at 
the Fermi energy projected on to each atomic sphere 
as a function of the distance from the interface for 
Co2FeAl/MgO/Co 2 FeAl. As can be seen, low spin- 
polarized interfacial states only appear on the first atomic 
layer. High spin-polarization is quickly resumed in next 
a few atomic layers. For example, the spin-polarization 
is recovered to 75% in the second atomic layer and it 
becomes 100% (i.e., half-metal) of 0.8 nm thickness. 

The non-equilibrium Green's function-based collinear 
spin transport calculations can provide a good de- 
scription of MR properties of conventional MTJs 
and good explanation on many experimental observed 



phenomenon [15|, Il6l . Il9j . However, the non-collinear 
transport method must be considered into the "per- 
pendicular" MTJs. We, for the first time, use the 
state-of-art non-collinear coherent transport method [see 
Supplemental Materials] to study the spin trans- 
port of CFA/MgO/CFA (001) p-MTJs. For com- 
parison, the widely-used CoFe/MgO/CoFe (001) and 
FePt/MgO/FePt (001) junctions are also included in this 
discussion. The calculated conductance values are listed 

in Table 1. As can be seen, Thus, CFA should be, at 

least in principle, an eligible candidate to substitute the 
conventional FM materials as the electrode material of 
the MgO-based MTJs from the perspective of enhanced 
MR performance. 

In conclusion, we provide a phase diagram of inter- 
facial thermal stability of Co2FeAl/MgO as a guide for 
experimentalists. A perpendicular magnetic anisotropy 
is demonstrated in CFA/MgO, which not only provides 
a magnetically thermal stability, but also reduces the 
stray field. An electric-field-induced giant modification of 
the PMA (a great magnetoelcctric effect) combined with 
a low damping constant and small stray field indicates 



a potential low-power and high compatible spintronics 
device application of CFA/MgO structures. We, thus, 
model a CFA/MgO/CFA (001) magnetic tunnel junction 
and find a high spin-polarized interfacial states and quite 
good magnetoresistance properties (large MR ratio and 
AMR) from non-collinear spin transport calculations. It 
shows that Heusler (CFA) based perpendicular magnetic 
tunnel junctions hold the promise for next-generation 
spintronics applications, such as non-volatile computer 
memories and logic chips [23] ■ We should point out that 
this prediction is just a proof-of-principle study, further 
experimental studies are needed. 
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